NPS ARCHIVE 
1968 
MCMILLAN, J. 


TEMPERATURE FLUCTUATICNS 
NEAR THE AIR-SEA INTERFACE 


by 


| John Gabriel McMillan 


LIBRARY 


NAVAL POSTGRADUATE SGHODS 
MONTEREY, CALIF. 93940 


UNITED STATES 
NAVAL POSTGRADUATE SCHOOL 


THESIS 


TEMPERATURE, FLUCTUATIONS 


NEAR THE ALR=SEA INTERFACE 


by 


John Gabriel McMillan 


June 1968 


TR Citatiori Database Page l-of2, . 


Add to Shopping Cart | Save for Later | 


Technical Reports Collection 
Citation Format: Full Citation (1F) 


Accession Number : 
AD0855533 
Citation Status: 
Active 
Citation Classification: 
Unclassified 
SBI Site Holding Symbol: 
AWS/TECH 
Field(s) & Group(s): 
080300 - PHYSICAL AND DYNAMIC OCEANOGRAPHY 
140200 - TEST FACILITIES, EQUIPMENT AND METHODS 
Corporate Author: 
NAVAL POSTGRADUATE SCHOOL MONTEREY CA 
Unclassified Title: 
Temperature Fluctuations Near the Air-Sea Interface. 
Title Classification: 
Unclassified 
Descriptive Note: 
Master's thesis, 
Personal Author(s): 
McMillan, John Gabriel 
Report Date: 
01 Jun 1968 
Media Count: 
53 Page(s) 
Cost: 
$7.00 
Report Classification: 
Unclassified 
Descriptors: 
(*SEA WATER, (*METEOROLOGICAL INSTRUMENTS, , SURFACE TEMPERATURE), 
DESIGN), THERMISTORS, CALIBRATION, INTERFACES, WIND, VELOCITY, POWER 
SPECTRA, PERIODIC VARIATIONS, THESES 
Identifiers: 
AIR WATER INTERACTIONS, *AIR SEA INTERFACES. 
Abstract: 
The design, construction and calibration of suitable instrumentation for measuring temperature 
fluctuations within a few centimeters of the sea surface, surface wave variations, and wind 
velocity are discussed. Small bead-in-glass thermistors mounted at 2.5 and 5 centimeters above 
the water surface were used to measure temperature variations under varying wind conditions and 
spectral densities computed on a digital computer. Spectra, co-spectra and coherence squares from 
0 to 5 Hz are presented. (Author) 
Abstract Classification: 
Unclassified 
Distribution Limitation(s): 


https://drols.dtic.mil/search97cgi/docview.dll?Key=AD0855533 &Query=%28%28memillan... 4/9/2003 


TR Citation Database Page 2-of 2+ » 


01 - APPROVED FOR PUBLIC RELEASE 
Source Code: 

251450 
Document Location: 

DTIC 
Change Authority: 

ST-A PER USNPS LTR, 1 MAR 72 


https://drols.dtic.mil/search97cgi/docview.dll?Key=AD0855533 &Query=%28%28memillan... 4/9/2003 


LIBRARY 


NAVAL POSTGRADUATE ScHoOl 
MONTEREY, CALIF. 93940 


TEMPERATURE FLUCTUATIONS 


NEAR THE AIR-SEA INTERFACE 


by 


John Gabriel McMillan 
Lieutenant Commander, Ufhited States Navy 
B. S., Naval Academy, 1958 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SClENCE IN OCEANOGRAPHY 
Erom the 
NAVAL POSTERADUATE “SCHOOL 


June 1968 


RN) eet) Ae pe ee ' 
M EDES 
oe, 
-T C./ 
OMT LLAN,S- 
ABSTRACT 
The design, construction and calibration of suitable 
instrumentation for measuring temperature fluctuations 
within a few centimeters of the sea surface, surface wave 
variations, and wind velocity are discussed. Small bead- 
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above the water surface were used to measure temperature 
variations under varying wind conditions and spectral 
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l. Introduction 

The energy exchanges between the atmosphere and the 
oceans are the driving force for a thermodynamic and mech- 
anical system which encompasses the fluid portions of the 
earth. The atmosphere derives energy from the sea in the 
horn Of sensible and latent heat and:in, turn, provides 
mechanical energy for oceanic circulations... Small-scale 
processes, such as convection and local wind fields, are 
also direct results of these energy exchanges, and the 
ability to accurately measure energy fluxes at any level 
over the sea would permit one to predict the growth and 
decay of such processes. 

There are generally four forms of air-sea energy 
exchange that one may consider: 

I. Electromagnerie radiabion 

2. Convection and conduction of sensible heat 


S. <CoOnvmeorlon andy diilttistom wr plerteme Ines) deme 
form of water vapor 


4. Vertical exchange of mechanical energy. 
In investigating the flux of any of these forms of energy, 
One MUSE ImeasUre Ginecthy certaim quantities. Por meat 
flux, the appropriate equations are: i =P C Tw and aaa 
PEW, whe re ES Lsvthe* sensible heat- flux; Ey is the latent 
heat flux, anda, are the densities of air and water re- 
spectively, © is the specific heat of air at constant 
pressure, L is the latent heat of evaporation (or conden- 


sation), T is the temperature and w is the vertical com- 


ponent of the wind velocity. Then, to measure the heat 


oe 


fluxes, one must be able to measure temperature, vapor 
pressure, and wind velocity at any desired level. Since 
waves on the water surface interact with the wind field 
(which carries the scalar temperature field), the sur- 
face wave field may also have an effect on heat flux, 

As a preliminary step in the determination of heat 
flux, methods of measuring variations of temperature, 
wind velocity, and water surface displacement on the 
smallest possible time and space scales must be devised 
and tested. This study was made to investigate some of 
the possible approaches to these measurements close to 
the air-sea interface and had as its objectives: 

1. The design of a suitable sensor system for 

measuring the temperature fluctuations within 
a few centimeters of the sea surface, the wave 
field at the same point, and the downstream 
wind velocity at a suitable level above the 
surface, 

2. The continuous recording of these three meas- 
urements in a form suitable for digital anal- 
VesnGs. 

3. The analysis of the records to determine the 
appropriate spectra and the relation, if any, 
between the air-temperature variations and the 
water-surface height under varying wind con- 


dit eELonms. 
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2. Design of the Measurement System 

A system to measure small temperature variations 
requires a temperature sensor that is durable, inexpen- 
Sive, and that has a sufficiently small response time to 
Eeecord Eluctuations at Erequencies of about 10 Hz. A 
small bead thermistor was used and although thermistor 
time constants are generally large compared to, say, a 
platinum resistance thermometer, it appeared that they 
were sufficiently small for this investigation. A Vic- 
tory Engineering Corporation bead=-in-glass thermistor, 
model 32A49, was used. The manufacturer's specifica- 
Hous fom ehismtimis are: a, nominal diameter of ©,013% 
Beime constant Cy) of 1.0 seconds im still air, a pow- 
SaEcdiccipacron Constant of 0.10 mw/ C, and a continuous 
operational range of -50°C to 300°C. The zero-power 


o 25%, and the resistance- 


Resistance at 22 CG is 2500 
temperature ratios are 5.6 for Re at CSG at 50°C, and 
ko) a a 25°C/R, at 125°C, However, over the range 
of temperatures encountered in this study (about 5°C), the 
resistance-temperature curve is essentially linear (Figure 
ioe 

Two of these thermistors were mounted on a guided 
styrofoam float, at about 2.5 and 5 cm above the water 
Surbace, as Shown in Figure 2. A variety of float ‘designs 
were tested for response to surface waves. The one used 


in the field was a hemisphere of about 2.5 cm diameter. 


In wave tank tests, movies showed the float to follow the 


waves quite well, although, of course, it tended to act 
as a filter for most waves smaller than the diameter of 
the float. 

To measure the water-surface fluctuations, a resist- 
ance wave gauge was mounted as close as possible to the 
thermistors, This gauge consists of a pair of 0.0L” @i- 
ameter stainless-steel wires, about 25 cm in length and 
positioned some 3 cm apart (Figure 1). A 5 volt, 20 Hz 
power source was used to avoid electrolysis, and the sig- 
nal output from the bridge was fed through a full-wave 
rectifier. 

Wind velocity was measured with a Casella sensi- 
tive cup anemometer, mounted approximately one-half meter 
above the water. The anemometer has a digital counter to 
record wind velocity and was used to determine mean wind 
speeds over time intervals of 10 to 30 seconds. 

The wave gauge and thermistors were installed as 
one leg of separate Wheatstone bridges and the outputs 
recorded on Mosley /7100B strip recorders. This recorder 
has individual DC amplifiers for each channel with maxi- 
mum amplification providing a 5-mv per 10-inch scale. 

The recorders were operated at'0.5 inches-per-second @hart 
speed for all measurements. The bridge and rectifier cir- 
cuitry for each instrument is shown in Figures 3 and 4, 
and a block diagram of the basic data-collection system 


Ls given in Figure 5. 


124 


The platform for field mounting the sensors is 
shown installed for measurements in Figure 6. The trans- 
mission lines for signal output are shielded, two-core 


telemetry cables with low-loss characteristics. 


3. Instrument Calibration and Observation Techniques 
The thermistors were calibrated by varying the tem- 
perature of a water bath and plotting the resulting re- 
sistance change for each temperature increment. The re- 
sults are presented in Figure 8. Since this curve was 
essentially the same for each individual thermistor, it 
was assumed that minor differences in slope could be ig- 
mored. Before field measurementsy the bridges were bal- 
amced im the Laboratory at a temperature close to that 
expected in the field. Each recorder channel was cali- 
brated to a scale of 1°C per inch and for equal sensi- 
tivity in each channel. The entire system was then test- 
ed with a known temperature change with the telemetry 
cables installed. The response time of the system was 
estimated by running the recorder at its maximum chart 
speed (2 inches per second) and recording a known temper- 
ature change. The time for 63% of the change (one time 
constant) was then measured directly from the record. 
Since the temperature changes were induced by adding 
water of a different temperature to the bathy the actual 
temperature change was not instantaneous, and the esti- 
mate of 80 msec is perhaps too large. It is, however, a 


goed andicabtion that “a flat response to 10 Hz-is a con-= 


| 
we) 
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servative estimate. This response will be somewhat less 
in moving air, and a response of 5 Hz was considered to 
be the upper limit for a flat response for this study. 

The wave gauge was calibrated in a wave tank to find 
the required immersion for maximum signal amplitude. The 
smaller the immersion, the greater the signal amplitude. 
Above a wave amplitude of about 5 cm, the output became 
slightly distorted due to the non-linearity of the gauge 
system. Again, the range of measurements in this study 
was sufficiently small that the response could be con- 
sidered linear. The output of the full-wave rectifier 
was compared with a common generated signal of O to 10 Hz 
and showed an undistorted output over this range. 

The Casella anemometer was not calibrated. The manu- 
facturer's conversion curves for counter readings to veloc- 
ities were assumed to be accurate, 

The field observations were made during the months 
of April and May, 1968 at Roberts Lake, Seaside, Califor- 
nia. This small lake is approximately one-half mile in- 
land from Monterey Bay and is about 450 yards long and 
200 yards wide. Between the bay and the lake are a series 
of sand dunes about 25 to 30 feet high which end about 150 
yards from the lake. There is no other obstruction to air 
flow in this area to seaward and the winds are essentially 
those present over Monterey Bay. The available unrestrict- 
ed fetch at the measurement site is between 75 and 400 


yards, depending on the wind direction. Most of the data 


was taken with a fetch of from 100 to 300 yards. The 
platform was placed in about three feet of water, 25 to 
30 feet from shore, and the signals were transmitted to 
the recording equipment with the low-loss cables. 
Measurements were made on six occasions, the first 
three with only a single thermistor mounted at 2.5 ecm. 
The lengths of the recordings were between 3 and 5 min- 
utes, with mean measured winds from 2.4 to 8 meters per 
second. The significant wave heights were between 2 
aad Chi, WLen periods Of from 172 to 12S "seconds “and 
lengths of from 10 to 15 cm. One representative record 
of the single-probe measurements and one of the dual= 
probe measurements were analyzed. A portion of a compu- 


ter drawn record of dual-probe data is shown in Figure 8. 


ay, Analysis 


In previous air-temperature measurements (Pond et 
al, 1966), a substantial portion of the spectral "energy" 
was concentrated in the region below 10 Hz. Since 5 Hz 
was considered to be near the upper limit of the flat 
response of the thermistor, the analysis of the data was 
restricted to the range of frequencies from 0 to 5 Hz. 

As the data was recorded in a graphic form, the conver= 
sion to digital form required a manual interpretation at 
imtervals of ©.05'' and ©.025°C.29This! preocesst@certaimly 
introduced some noise into the record (in addition to 

that generated by the electronics in the system), but this 


was unavoidable. The ensuing analysis showed this noise 


igor 


to be negligible. 

Since the temperature fluctuations at the air-sea 
interface are felt to be essentially random and station- 
ary, it is appropriate to analyze the data by computing 
power spectra. This method of analysis is the basis of 
an existing FORTRAN IV program in the library of the 
N.P.G.S. computer library, PROGRAM BLACKY (Table 1). 
This program computes smoothed power spectral estimates 
from the Fourier transform of the covariance for a given 
frequency bandwidth. 

In analyzing any continuous process using discrete 
points, the problem of aliasing is encountered. For this 
study, it was assumed that little energy existed above 5 
Hz. PROGRAM BLACKY also proviced a filter by running a 
weighted average for three consecutive points of the 
record. This would reduce the possibility of folding and 
also remove most of the high frequency noise. 

In order to confirm the results of PROGRAM BLACKY, 
and to obtain a graphic display of the spectra, the data 
was also run in PROGRAM BIMED O2T, an autocovariance and 
power spectrum program of the Health Science Facility, 
UCLA. The records were also analyzed for means, density 


and cumulative distributions. 


2. Results and Conclusions 
The data was processed for smcothed estimates of 
spectral) density. froma0 Go Sj jhzyat ansimeerva LBoLAeo02s 


cycles. The record was filtered by means of a Hanning 


— 
ion) 
eo 


lag window of the form: 


where De is the weighting function, m the number of lags 
Gor the secord, and © is an imcremental index of m.. The 
weights for this analysis were 0.25, 0.50 for the three 

consecutive points. The first and last points were com- 


puted with weights of 0.5 as (ORS) x end (0.5 ox + O59 


m=1 
Xn This low-pass filter effectively removed the high 
frequency noise of the record. Since little energy was 
felt to exist at high frequencies, the loss of this por- 
tion of the signal was considered to be insignificant. 

The smoothed spectral estimates for the temperature 
variations are presented in Figures 9 and 10, and the 
corresponding wave amplitude spectra in Figure ll. A co-= 
spectrum for the dual-probe temperatures is presented in 
Figure 12. The values from which these plots were taken 
are given in Tables II and III. The spectra are presented 
as log, o¥(£) versus logj(£); Wa tle SAG@iD ein) cians) sons 


(°c)*sec. The spectrum is defined such that 
ao — 
fv (eae = 6 
° 


where 9 is the temperature variation (Batchelor, 1959). 
The wave amplitude spectrum, @(f), is defined similarly, 
with units of cm“sec. Since Taylor's hypothesis of “fro- 
Zemectucpubence (laylor, 1938), is net valid fer a recion 


of shear flow, the relation 
Zeit 


where k is the wave number, f the frequency, and U the 
mean wind velocity, was not used and the spectra are pre- 
sented in terms of frequency rather than wave number. 

The total "energy" was estimated in two ways in PRO- 
GRAM BIMED 02T; (1) by computing the variance, and (2) by 
integrating the power spectrum with the trapezoidal rule, 
The ratio of (1)/(2) gives the fraction of the total ener- 
gy in the computed spectra. The following table gives the 


individual values and the resulting ratios: 


ae ee eee ere ee ee ere errr ore ee ee 


Y(£)/B(£) Check sum Energy by Brees Ratio 
ae seo Jijuei Of 2PSDs _poTrep. ule av {aisle SIGTIOCeee 
Temp. 2.5 cm 1.1057844 1.1058750 -9.06x10"? 1.00 
(Single-probe 

Temp. 2.5 cm 0.2300338 0.2300585 2.47x107> 0.99 
(Dual-probe ) 

Temp. 5 cm 0,3806309 0.3806635 3.27x107? 0.99 


(Dual-probe ) 


Wave Amplitude 2.5749149 2.5751486 2.337x107°* 0.99 
(Sing le-probe ) 


Wave Amplitude 0.29007 39 « 2.2900738 100x107? 1,00 
(Dual-probe ) 


The high percentage of energy in the computed spectra indi- 
cates that there is little energy above 5 Hz and there was 
negligible energy folding from higher frequencies. 
Confidence limits for the spectral estimates were 
calculated for an 80% interval, For the analyzed record 
length of 120° seconds, a sampling interval of 0.1 seconds, 


and a correlation.lag value of 200, the 80% limits are; 


0.65 W(f) ¢ W(£) ¥ 1.91 (£). The computation of these 
limits is presented in Appendix A. 

Coupling, between the surface waves and the tem- 
perature £Lluctuatbions at fixed heights above’ the undu- 
lating surface, does not appear to exist, even at the 


higher wind speeds. The coherence square function is 
2 


2 Cea Ge) 
y (f) = “1 ‘ 
XY Cece) G,(£) 


where Cu) is the co-spectral density function and Gag) 
and G,(£) are the individual spectral density functions 
(Bendat and Piersol, 1966). The coherence square for the 
temperature and wave spectra (Table II) shows virtually no 
Correlation between Ehese quantities.  Simee Ehe tempera- 
ture measurements were taken at a height that remained 
Constant with respect to the water surface, this lack of 
correlation suggests that the temperature structure at 
this level is not affected by the passing wave surface. 

In contrast, spectral densities for temperature measure- 
ments taken at a fixed height above mean sea level have 
shown periodic maxima at the same frequencies as wave 
amplitude spectral maxima (Makova, 1963). A comparison 

of measurements taken at a fixed height above the surface 
and simultaneously at a fixed height above mean sea level 
would be useful in further studies of this problem. 


The one-dimensional spectra of temperature fluctua- 


Ee 


PLoms ported on a Loe-loo scale (Bieurces 9 and 10) show 
Pp & & & 


a =e slope for frequencies above 0.2 Hz. This slope has 


(et 


also been found for temperature measurements taken with a 
platinum resistance thermometer at heights of one meter 
and greater (Pond et al, 1966). Since the inequality, 
kz 277 4.5, does not hold near: 0.2:Hz, ‘this result indicates 


that asrelation of the form 


Ly 5/ 
PCE Ste hee Stk Fg 


8 
(Batchelor, 1959) is valid for anomalously large scales. 
(Here, GC) is a conversion constant to relate the spectrum 
to frequency rather than wave number, eg is the scalar 
dissipation, e is the energy dissipation, and A is a uni- 
versal constant.) This result was also found by Pond et 
al (1966). The high values of the water-wave spectra be- 
low 0.05 Hz are not readily explainable, but may be re- 
lated to a periodic pile-up of water at the leeward end 
of the fetch, where the measurements were taken. Ina 
small lake, such as was used in this study, periodic 
gusts can easily set up a seiche. 

An examination of the coherence square for the dual- 
probe temperatures shows a high degree of correlation over 
the entire spectrum (Table III). This, together with the 
similarity in the amplitudes of the fluctuations as seen 
in Figure 8, is somewhat surprising in view of the 5-cm 
probe being twice as far from the water as the lower 
probe. The mean temperatures for the dual-probe record 
aber a4 26°Gvaity 275 cm and 14 .04°C at > cme ihe mean 
air temperature as recorded at approximately 5 cm above 


the water with a mercury thermometer was 14.0°C and the 


DOr. 


water temperature at the surface was 20.8°C. Clearly, 
the air near the water was highly unstable, and it ap- 
pears that the greatest air temperature gradient occurs 
below 2.5 cm. 

The results of this study indicate that the data 
collection system used for obtaining temperatures and 
wave records was reasonably accurate and that additional 
work should be done to reduce the level at which the tem- 


peratures canbe taken. 


6. Recommendations for Future Studies 

Further work in this field should include a con- 
tinuous record of the wind field, additional thermistor 
probes at lower levels (and at least one below the water 
surface), a simultaneous record of temperature at a fixed 
height above mean water level, and the development of a 


sensor mounting that will not interfere with the wave 


field, such as a servo-follower system. 
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APPENDIX A 


Confidence Interval Calculations 
Record length = LS 120 seconds 
Sampling interval = f = 0.1 seconds 
No. of lags per record length = m = 200 


Degrees of freedom = d 


No. of data points = N 


eee (Blackman and Tukey, 
1958) 


iS x ry (Bendat and Piersol, 
Kqse/2 Xa; 1-%/2 1966) 


Assuming that the Fourier transform bears the same rela- 
thom, Of = 0.2 and d-= 12, eives values for 9 ener and 
Xa;1-4/2 a 

2 
Xa;%/2 = 26.22 
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